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SECTION I 



The adaptability of the differential analyser in solving 
the turbo prop control system problem was investigated in this study, 
and some of the more fundamental results which are obtainable by this 
method of approach aro herein presented. 

The basic control assumption used in this investigation was 
that rotational speed be governed by blade angle. Temperature was 
neglected, and fuel flovr was considered as o step input resulting in 
an instantaneous increase in torque. 

It was found that adequate response at a given operating 
condition could be obtained by essuning blade angle a function of off- 
sneed find integral cf off speed, i-.t other operating conditions, however , 
there was a wide variation in response. This variation in response was 
reduced by requiring blade angle to vary , in addition to the previous 
assumptions, as the derive tive of the off speed. 

It v/as concluded that the differential analyzer has definite 
merit for application to this control problem, both for analysis and for 
possible actual installation as the c ntrolling unit. 

This investigation is but an infant's atop toward a very adult 
problem. Additional and amplifying procedures and methods of approach, 



not possible of accomplishment due mainly to lack of time, are 
hei'einafter presented. 



SECTION II 
nrmoDUCTiQN 



Analysis of the characteristics of oov/er plants available 
for the propulsion of aircraft clearly indicates the place for the 
turbo-prop paver plant. The highly developed reciprocating engine 
has proven itself in tho speed range up to 400 MH 0 The turbo-jet 
would seem to be the answer where very high speed, 550 I xH end above, 
is specified; but this high speed is accompanied with the disadvantage 
of short range due to the high fuel consumption of the turbo-jet engine. 
The turbo-prop power plant fills the speed gap between 400 end 550 i.lHi 
and in addit ion, excellent range can be obtained, ’.Vith the development 
of the supersonic propeller, the speed of turbo-prop powered aircraft 
might well reach that of trrbo-jet powered craft. Thus the need for 
the turbo-prop power plant is apparent, both for commercial end mili- 
tary operations. 

The turbo-prop power plant consists essentially of a gas 
turbine driving a propeller through reduction gearing. The control 
problem for such a system is inherently complex because of the number 
cf variables which enter into the problem end the interrelation of the 
variables. The five primary variables are Engine Speed, Torque, 



Temperature , Fuel Flcnv , and Propeller Blade Angle. In addition, 
maximum engine efficiency is attained at nearly maximum allowable 
turbine speed and temperature, thereby requiring operation in a 
regime havirg narrow oontrol limitations on certain of the variables. 

The ultimate requirement for a turbo-prop control system 
is that the pilot be able to select a desired value of power, without 
exceeding any of the design limitations of the power plant components. 

A great amount of research is at present being conducted in an attempt 
to obtain such a control system. 

In References 1, 2 and 3, a general discussion of the control 
requirements for turbo-prop power plants is given. In Ref. 4, the 
adaptability of the electronic differential analyser to the analysis 
of the requirements of the control systan is discussed. In Ref. 5, the 
utility of the electronic differential analyzer for obtaining solutions 
to engineering problems is considered, and a complete description of the 
type analyzer used in this investigation is given. 

It was the purpose of this investigation, in addition to giving 
the authors practice in the use of the electronic analyzer, to carry out 
a study of the basic requirements of the control system for a turbo-prop 
power plant. Since the purpose did not include the design of such a 
system, a fundamental approach has been attempted throughout. The method 
follcr./ed was to start witti an initial, simplified equation representing 
the basic characteristics of the turbo-prop power plant, and to add to the 
initial system, as time allowed, certain functions which would improve 
the response characteristics of the control system. 

It is regretted that due to time limitations, a more complete 
investigation could not be raadeo 



The authors wish to express their appreciation to Doctor 
L. L. Rauch, Department of Aeronautical Engineering, University of 
Michigan, for his assistance and guidance throughout this investi- 
gation; and also to Professor J. //. Luecht, Department of Aeronauti- 
cal Engineering, University of Michigan, for his assistance in 
obtaining the engine and propeller data used, and in the reduction of 
t;ds data to a usable form 0 
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SECTION III A 



DEYELOH&IIT OF BASIC STATIONS 



Definition of symbols end terminology: 



Subscripts: 

a - aerodynamic 
c - compressor 
e - engine 
p - propeller 
t - turbine 
\{f- fuel flow 

Superscripts: 



~ ( ^ Oc y ^ ^i)jL. e 

C fji = 

C f»At = 



I - Moment of inertia of the engine 

G 



- steady state 



Symbols: 



C, 





to the (engine 



Q, - Instantaneous aerodynamic torque absorbed by 
& 

the propeller 

Q - Instantaneous net engine torque (equals CL-Q. ) 

Q^ e - Instantaneous engine torque delivered to the 
reduction gears 

Qp - Instantaneous torque to propeller from reduction 

gears 

( ■* e “-<p e ) “ Instantaneous torque available for acceleration 
of the engine 

(0 -Q. ) - Instantaneous torque available for acceleration 

rr 0. 

of the propeller 

q => Qp e “-^ e “ Difference between instantaneous and steady 
state torque to reduction gears 
i = - Reduction gear ratio 

JX - Rotation speed 

J - Blade angle 

9 - New variable Introduced for ease of mathematical 

treatment 



Relati ons: 

Transient 

% • V 1 

«e ’ - % 

%-%e- X e 

V 4 a * : p (d V dt) 

Also; d-rt^/dt = d(-/le -yig)/dt 
dn^/dt = &{-Ap -Jip )/dt 



Steady State 



* 

>e 


V 1 






= 




% 6 


q* 

pe 






a 


TPP 



since steady state 
derivative is zero 
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Development of engine equation? 

Ilote: Constant density ratio, temperature ratio, 

and flight velocity assumed here and also in 
the development of the propeller equation » 

From above : 

Qe-Qpc. = x e c *( yi e ) 

Q e (v4,^-e)= +dJLe) } 

0* e + OQ*ft w f)j t£ ^ Wj r + ^ yL£ ) 

Q*e +C + C ^ £ djLe - 

* * 

Q% + ) <Kno = ® ? e * 

Therefore, the engine equation is: 

x e JlJL'-jft/Jt = f * C tjLe -9. (1) 

Development of propeller equation: 

From above: 

Q n - <?o. = L p U(ji P --n*e/Jt . 

Bierefore ^ . 



where 

<?e = 
O* - 
O e = 

But 

<5(*d “ 



where, ^ = Q a (-*£,/) = ; 

Qj. ~ G«. ** (J<2<*/ J-Ag^ J-rt-C -f Ofly/ g 

<?•- = <?*. •+ -*• ** •* ^ j 

<?*>£ 3 + £/•** **■ 

Therefore , 

Zf>JjL x el (*e-JLe*)/Jt - Qfa ~ QIA ~ C fre t * Jle 

But » # ♦ * 

fye = tye +y ; a-** ' 

Therefore the propeller equation is: 

X ?/JL 31 ~ t ~ C fi*e Jjle ~ C & (2) 

Development of combined equation: / A \ 

Te <A( JL e-A*)/At - C *w f J y f + c ^ JLe " fr- 

If/Ji «t (Ae -JLe )/jit. ~ T ~ < ^ JLa - Cpp Afi , (£ ) 

Eliminating q and substituting differences for differentials: 

(J e * x r/jA) <l( JL e-JLe*J/lt + c fj3&y3*) 

+(Cf Ae +c< Jtt )tf-£) = c t " f (*A). (3) 

Equation (3) is the basic equation used in this investigation. 
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SECTION III B 



? 

DATA ON THE T. G. 100, TURBO-PROP ENGINE 



For purposes of this investigation, operating data on the 
General Electric T. G. 100 turbo-prop engine, as presented in Ref. 6, 
was used exclusively since this was 12ie only data available to the 
authors in an unrestricted form. The fact that the T. G. 100 is 
essentially a constant speed engine, does not reduce its value as 
representative of a turbo-prop engine, but does simplify the control 
probl an to some extent 0 In a variable speed engine, a new engine speed 
accompanies each new power plant output called for, in order to operate 
at maximum system efficiency. This affects the control system problem 
in that this new engine speed ?/ould have to be called for, or scheduled. 

For this investigation, a flight velocity of 300 MPH was 
assumed for Step I; and for Step II, velocities of 100, 300, and 500 
M EH were assumed; all at a fuel flow rate of 1400 lbs/hr and at sea 
level. In order to evaluate tiie coefficients for the equation repre- 
senting the power plant, the values of total shaft horsepower listed 
in Table I were obtained from Figure OA of Section 24, Ref. 6, repro- 
duced in this report as Figure 2. The friction loss was obtained from 
Figure 23 of Section 23, Ref. 6, and subtracted from the total shaft 



horsepower thus obtaining the Engine Shaft horsepower This 

was tion converted to engine torque ( ^ ) . 



[ Qe ' fL 



Lbs = SJ-R, * 



33, OOO 



] 



A plot of Engine Torque 



e " zn(n e) Rpr\) 

( s, ) versus Engine Speed (J"l e ) with fuel flow rate ( v/ ) as parameter 
wes then made for 100 , 300 and 500 L2FH flight velocities in Figures 3, 
4 and 5 respectively. The slope of these curves gives the coefficient 

for the operating condition (vel city, and eng ne 
power) chosen. The value of far flight velocities of 100, 300 

aid 500 !.TH at sea level and w^ of 1400 lbs/hr are given in Table III. 

The polar moment of inertia of the rotating elements of the 
engine (turbine, compressor rotor, reduction gearing, and accessories) 



r \ 

erL e v d -fl e 



is given in Hef. 6 as I g = 4.008 slug ft.‘ 



SECTION III C 



3ELECTI ON uF IEOPELLER 



Definition of symbols: 

P 

Cp (Power Coefficient) = 

^ 3 mas3 density of air (slugs per cubic foot). 

/7J. 3 propeller rotational speed (revolutions per second). 

0 = propeller diameter (feet). 

\/ = airspeed (feet per second). t 

- blade angle (degrees) at 0.75 propeller radius. 

Since there was no available information on the propeller 
intended for use with the T. 0. 100 engine, it was necessary to select 
a suitable one. The best information available on propellers suitable 
for this type installation /as found in Ref. 7. The theory as presented 
in Ref. 8 was utilized in conjunction with Ref. 7 to determine the most 
suitable propeller possible. 

Ref. 7 contains information only on ten foot diameter propel- 
lers. Of the propellers considered, the one finally determined as 
capable of absorbing the engine power available was a six wide blade, dual 
rotating tractor propeller, the characteristics of irhich are given in 
Figure 16 of Ref. 7, reproduced as Figrre 6 of this report. In order to 
check the suitability of this propeller, combinations of flight speed, 



engine power, and engine RFtl were selected at sea level, and — - and 

/l\ i/ 

C a computed. Entering Figure 6 with these values of ^L. and Cp , noints 

r /*./> r 

indicating j£ and propeller efficiency were determined. Plotting a 
number of these points for speeds of 100, 300 and 500 LUK indicated that 
at 300 HPH the points were to the left of maximum efficiency, While for 
500 1.PH, the points were to the right. This phenomena was caused by too 
large a variation in JC due to the fact that an initial basic premise 

/K. D 

was constant engine rotational velocity for 'steady- state operation 
This loss of propeller efficiency was disregarded because, for the pur- 
pose of this analysis , only adequate power absorption of the propeller 
is important, not power output. Since this situation does exist, however, 
it is implied that even though an engine of itself gives maximum economy 
at full rated RFi£, a simultaneous schedule of fuel flow and RFM instead 
of fuel flow only, is necessary to compensate for the loss in efficiency 
of the propeller. 

Having decided that the propeller was acceptable, the arguments 

V 

for Table II were obtained from Figure 6 by entering this Figure with^jj-p 
and reeding of f Cp for various . From this table, Figures 7 through 
Figzre 12 were constructed for the three speeds selected for analysis. 
These curves are plots of ^ vs B with -A e as parameter and VS -fl-e 
with JZ as parameter. The puipose of these curves is simply to enable 
determination of C ^ and for each operating condition analyzed so 

that these coefficients might be inserted in the general differential 
equation representing the system. and are slopes at the par- 

ticular operating condition, being equal to 

respectively. The coefficients for flight speeds of 100, 300, and 500 
HPH at sea level and w of 1400 lbs/lir are tabulated in Table III. 
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One other quantity needed is the moment of inertia of the 
propeller, .diether of the solid aluminum alloy or hollow steel type, 
the result was essentially the same. One method of approach is illus- 
trated in Ref. 6, and is applicable to aluminum alloy blades. The 

2 

estimated Ip using the method of Ref. 8 was 70.2 slug feet . 



\ 
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SECTION HID 



THE ELECTRONIC DIFFaHEIjTIAL HELriLR 



A complete, detailed description of the electronic differ- 
ential analyzer used in this investigation is given in Ref. 5, end 
therefore vail not be repeated herein. In general the analyzer con- 
sists of plug-in feedback a^lifier units vvhich by simple external 
changes in connection can be used for integration, differentiation, 
multiplication and division by constants, summation end sign changing. 

For recording data, the desired outputs from the analyzer 
./ere fed into Brush, J/odel BL-913 direct current amplifiers (modified), 
the -utout of which ;vas f--d intc a t. to- channel Brush, iiodel BL-20C, 
magnetic tyoe, recording oscillograph. 

A picture of a typical analyzer set-up is given in Figure Id. 



SECTION HIE 



PROCEDURE 



The basic equation representing tiie power plant is equation 
(3) (see Section III A above). 

fx e + C-^e ~-^t) t is 6 ?* ) + IW C ejl £ ] (-*■£ ) 

= Ce„ # ("f ' "S) . (3) 

For purposes of analysis, a new variable is introduced, i.e.; 

a = O 1 * ~JT-e) } 

e -- (** . 

Substituting © in equation (3), the basic equation becomes: 

= Ce„, (»'♦-"£>. (^j 

If it is now assumed that the propeller pitch control changes 
propeller pitch at a rate proportional to the engine off speed, thm 

*^4tt ~ h (-Ae-Sie ) } 

J = J>* + k t £(Ac~Jle)Jt =J*+kiQ } 

J&-JS* a k,J (-Ae-Jll) s kt & ) i'S) 



where is "tb.© governor proportionality constant. This type of 



control on propeller pitch will he termed "integral control on J2 
Substituting this value of (J$ -y£* ) in equation (4), the 
following second order differential equation is obtained: 



Figure 12, indicates that the coefficients of equation (6) remain prac- 
tically constant over a reasonable range of operation, the inertia term 
naturally remaining constant . 

If the coefficients are considered as constant, equation (6) 
is comparable to the standard single degree of freedom vibratory equation 



The known procedure for the analysis of this type equation can now be 
followed.- 



of the system is insuf fici ent for adequate response. 

The natural damping was therefore augmented by requiring the 




Analysis of the engine and propoller curves, Figure 3 through 





is the inertia or mass teim, 




damping terra, 



and 




is the spring constant term. 



(A) Step 1. 



As indicated in Ref. 3, the natural damping term 




prop pitch control to change pitch at a rate proportional to the 



derivative of engine offspeed, in addition to the engine offspeed 



J-A/Jt = k, (A' -jU) + ki, (-A6 --A?) , 

^ * * 

= A,® + k A & . 



The basic equation after substitution then becomes 

[z* * wj « + [ C tA< - <*« + % k*] ® +fy k ‘ e 

= C' Wf (*i-vh. 

For purposes of analysis, it is convenient to consider the 
right side of equation (8) as a step input of engine torque (q^ ) 
resulting from a step input of fuel flow rat® ( ; and then to 

determine the response characteristics of the system due to such a step 
input. It is considered that this is a conservative approach, since 
the effect of an actual change in fuel flow rate, involving necessarily 
a finite rate of change, would represent a less stringent condition. 

The basic equation, equation (8), can now be expressed as 

[*« ' T r/J} ® - C ejl£ + C fji t,] e + C e j k, e 

~ - ) 

whore ~ fi.+ VW is the inertia term, 



a l, 




is the augmented damping term. 






[ A, J 



is the spring constant term. 
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Various values of damping rutio { ■ ) were next 

assumed; and for each value of damping ratio, various values of the 
natural undamped frequency ) of the system were assumed. 

Since the value of the inertia term (a ) is known and remains constant, 
the value of the damping term (a-^ ) and the spring constant term (a Q ) 
could be determined for each combination of ^ and <4^ chosen: i.e, 

^ , I e Ip/+ A = H-ooV +.SH X - } 



CL, 






'/*. 



= V. SS <4* 



a~ l = 4. ^ • 



For this part of the investigation, a cruising operating 
condition at sea level was assumed, with a flight velocity (V ) of 300 
IJFH and fuel flow rate (w^) of 1400 lbs/hr. For this .condition, the 
values of C *l£ , c e Af t and Z 1A had previously been determined and the 
values listed in Table III. It was then possible to evaluate the pro- 
portionality constants k^ and kg, since 

k > 3 a °/ c jv , 



Aj. = -hi- ~ ' C f** 



- c 



ejie 



't* 



Far completeness of information, the value of tho natural 
damred frequency (fa)) of the system for each combination of and <4* was 
determined by u) = u)^ J / - . 

The values of a . a., k, , and k 0 for the combinations of ^ and 
chosen are given in Table IV. Tho reciprocals of the coefficients 
are also given since they ore used in the analyzer set up. 



All coefficients having been determined, the differential 
analyser was set up. The diagram of the analyser circuit is given in 
Figure 18. Amplifiers A^ through are used to solve the basic 
equation (equation (9j; while amplifiers A^ and Ag were added to solve 
the engine equation (equation ( 1 )) for q, in order to record this variable. 
The step input of torque due to change in fuel flow rate ) was applied 

as a battery voltage across the inputs to amplifiers Ag A 7 when Ihe 

switch S was closed. An example of the record of a typical step input 

) on the Brush recorder is given in Figure 19. 

In order to stay within the voltage limitations of the ampli- 
fiers (approximately 0.5 volts to 100 volts), it v/as necessary to change 
certain input and feedback amplifier resistances by factors of ten. 

This resulted in some of the variables recorded being different from the 
actual values by a multiplication factor, the factor being indicated in 

Figure 18. A correction for this factor was made when the recorded data 

was reduced. 

ic 0 

The outputs of amplifiers A^, (Fi-£--h-e ); (0); Ag (^-pS ); 

and Ag, (q) were recorded. 

For each combination of $ and *< 4 *, the proper coefficients, as 
given in Table IV, were inserted in the analyzer in the proper position, 
as indicated in Figure 18. 

A step input, , v/as then applied for each combination and 

if 

Brush oscillograph traces taken of the response of the variables (-Ae-Jl 6 ), 
(q), iyS-yQ*), and (©). The recorded oscillograph traces are shown in 
Figure 14 through Figure 17, respectively. 



(B) Step 2. 



As will be 3een from the results of this step, the response 
of the control system varies over a -wide range with changes in flight 
velocity if the type of propeller p£tch control described in Step 1 is 
employed. If for instance, the proportionality constants k^ and k g 
are chosen to give proper response at 300 j.PH , then these same constants 
coupled with the coefficients corresponding to a different speed cause 
less damping at lower fligit velocities and more damping at hi^ier 
velocities. This change in response with change in speed is illustrated 
by the 0$ D.C. curve of Figure 30. (# D.C. refers to the percent deriv- 

ative control used and is defined later on in this section). It would 
be desirable if the control system could provide approximately constant 
optimum response over the complete range of flight velocity. To 
approach this result we will therefore require the propeller pitch 
control to change pitch at a rate proportional to the second derivative 

of the engine offspeed in addition to the requirements of Step 1. Thus 

- k, (-A-e ~-Ae) + kj,(-A-e -JL £ ) + (-A-e -ji e ) ^ 

+ k x (A £ -JLe*) + Aj (-A* '-rt-e) } 

- k, Q + k±e + ki e f 

where k„ is another propeller pitch control proportionality oonstant. 
o 

This added term can be called "anticipatory" or "derivative control on 

JB -• 

l 
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Again substituting in equation (5), the basic equation for Step 2 



becomes 

[z« + %•> + c fji *»} a + i c ^ c ■* Cfrkje 



where 



+ Cfyk, Q =» 

*a 3 L 1 * + Ip /i x + C K* k s) 
£L ‘ “ L C i^B + ^ 

< 2-0 - 7 




is the inertia term, 
is the damping term, 
is the spring constant term. 



It is seen that the addition of derivative control on Jt can 

i 

be considered as an addition to the inertia term (ag) of equation (10), 
thus giving a new "effective" inertia. Hie basic equation remains in 
the form of the standard vibratory equation and can be analyzed in a 
manner similar to that used in Step 1, 

An undamped natural frequency (<*^) of 5 rad/sec was now arbi- 
trarily chosen as representative of this type power plant at an operating 
condi tion of 300 LTH, and a damping ratio of 0.6 was chosen as an 
optimum overall value at this speed. Since the derivative control on 
Jb can be considered to change the inertia term (ag) of equation (10), 
the procedure followed was to increase the value of the inertia term 
by convenient percentages, i.e. 100$, 200$, and 300$. This percent 
increase was termed percent derivative control, and is abbreviated here- 
after as "$ D.C. " 

In order to cover the complete range of flight velocities, 
three representative operating conditions wore chosen, i.e. low velocity 
of 100 MPH, cruising velocity of 300 and high velocity of 500 MPH; 



ail at sea level and a fuel flow rate of 1400 lbs/hr u The determina- 



tion of the values of the coefficients , and for these 

operating conditions has already been described, und the values ar 4 * 
given in Table III. 

By maintaining ^ *= .6 and <0^ » 5 for all values of $ D. 0. for 
the 300 J.IFH operating condition, it was possible to determine the values 
of k^t kg, k 3 , and & 2 , a 1 , a Q for each case. 

[ To oc][ic +_ *>A*7 

k3 c 

a 2 = [/ + %0c] L?c + T ^ x ] 

a 0 3 °-X 



k 



1 



<X -°/ 



a x = A ^ 7^7 = £ C />_A £ ~ + ^ 




f>JL£ + ^e-^-6 

C fjZ 



Hie calculated values of a^, a^, a 2 , k^, kg, and kg for flight 
velocity of 300 MEH are given in Table V. 

The values of k^, k 2 
for the 300 1TH operating condition ware used for the corresponding % D.C. 
at the other flight velocities chosen. It was than possible to calculate 
the coefficients (a o , a^, Qg) of the basic equation at each % D.C. and 
airspeed, and therefore the ^ and u) . 



, and kg thus determined for each % D.C. 



a o “ C fJ kt 

3 C f-n-e ~ Ceji £ + 
a „-/"/-#■ %/>C ] [ X* 1- T f>/^] 

Ci ^ 

5 » ^ / / / * /«-o 

The values determined for these two operating conditions 
are also given in Table V. 

In Figure 20, a plot was made of $ 73 flight velocity (V Q ) 
to show the theoretical effect of derivative control on the damping 
rat io . 

The coefficients of the basic equation having been determined 
for all caabinations of operating conditions and % D.C. chosen, the 
response of the system to a step input of torque (q^) was obtained using 
the malyzer in the same manner as described for Step 1, the some vari- 
ables being recorded. Hie modified analyzer circuit diagram is given in 
Figure 18. The recorded oscillograph traces ere shown in Figure 21 
through Figure 24. 
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SECTION IV 



RESULTS ,ND DISCUSSION 



From the oscillograph of the four variables recorded, the 
maximum aid steady state deviations in millimeters, slope in. the case 



ments in millimeters were translated into voltages by applying the 
calibration data and attenuator settings on the recorders, and also 
the multiplication factors made necessary by the voltage limitations 
of the analyzer amplifiers. These voltage magnitudes wero then equal 
quantitatively to the value of the variable measured. The values of 
frequency and time to damp were difficult to obtain from the traces 
and tlierefore the absolute values given can be considered only approxi- 
mate. Values of frequency not listed in the tables were impossible to 



measure from the traces. An example of the method used in the reduction 
of data is given in Section VII, Sample Calculations. 



65$ power® The step input of torque was 63.5 ft. lbs. or equivalent to 
roughly a 10$ change in power based on the original steady state power. 
This corresponds to about a 6.5$ change basod on maximum power. Since 




time to damp were obtained® These measure- 



Step 1. Results 



Step 1 was based on an original steady state torque of about 
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the coefficients of the differential equation for the system change 
only slightly with change in output of the engine, the equation can 
be assumed linear for all values of torque at one airspeed without 
too much error. For instance, if a jump from 65% to 100% torque v/ere 
desirable, the results tabulated oould be multiplied by 100-65/6.5*5.4 
to give fairly good approximations of response. 

In this step, it was assumed that^S was a function of the 
integral of the offspeed and the offspeed. The traces taken are pre- 
sented in Figure 14 through Figure 17, and are arranged to show the 
effect on response of changes in * and co . Only relative and not 
quantitative comparisons of maximum and steady state values of and 
9 shovm on the traces may be made, due to the necessity of having dif- 
ferent attenuation and calibration values set on the recorders. This 
was necessary in order to get traces of reasonable size. The reduced 
data is given in Table IV, and Figure 25 was drawn from the tabulated 
values. 



Figure 25 shows the effect on maximum response values of q, 
and ^ of varying damping ratio (^J) and natural undamped fre- 
quency It is seen that increasing results in lower maximum 

values of torque and engine offspeed, end higher maximum values of ji 

and 4^ . Increasing , in general, decreases the maximum values of 

It J 

the variables plotted except for in which cese the reverse is true. 

at 



In an actual control system, there would naturally be limits 
on the maximum allowable or obtainable end • These limitations 



would partially define the design problem which oould be solved by 
selecting proportionality constants (k^ and kg) to obtain the necessary 
^ and t * J n for adequate response. 
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wa s not considered previously in this investigation 



Jt* 

dt± 

due to the inability of measuring this quantity in all cases from the 
A traces. Several approximate values were obtained from the ^ = .3 
traces; 




These average values 'were obtained from high speed traces 
(not included) by dividing ^^') t1AIL ihe tine interval necessary to 

J 

All values of could have been obtained in this 



reach 



(*J*\ 



fashion using the differential analyzer by merely changing the time base 
of the system or by setting up the system to record this varirble. 

It is understood that a typical system now in use would at 



least be capable of a 






of 10 /sec. in .03 sec. and therefore a 



A*a o , 2 

of 333 /sec . Such values, nee determined fren the propeller 



governor ctuoracteristics , could be set on the differential analyzer. 



Step 2. Results 

Step 2 was based on an original steady state torque of about 
65;o power at 300 j*PH. The step input of torque was 62 ft. lbs. or 
equivalent to roughly an 8 ,* change in power based on the original steady 
state power at 300 KPH. These percentages varied slightly with change 



in speed, due to maintaining constant fuel flow instead of constant 
torque throughout the speed range. The comments made for Step 1 
regarding linear assumptions apply here also. 

In this step, derivative control on^ wa3 added to the pre- 
vious requirements of the propeller pitch control to compensate for 
changes in systm response with changes in flight velocity. 

Prom Figure 20, it may be seen that, theoretically, damping 
ratio, end therefore system response, can be made to approach a con- 
stant optimum value for the complete range of flight velocities by the 
addition of the proper mount of derivative control. The oscillograph 
traces for this step are presented in Figure 21 through Figure 24, in 
a form which shows the effect of changing flight velocity at any c /o D.C.; 
and also the effect of increasing the D.C. at any chosen flight 
velocity. It is seen from these figures that the rosults agree with 
those predicted theoretically. 

The reduced data for this steo is given in Table V, end 
Figure 26 was drawn from the tabulated values , 

Fran Figure 26, it is seen that by increasing the £ D.C., the 
peak displacement cf +he virialles is decreased for all flight velocities. 
The advantage of 3uch attenuation are apparent. For instance, by decreas- 
ing the transient overshoot of torque, excessive reduction gear loads can 
be prevented. 

It is noted from the response of . /3 , as seen in Figure 23, 
that very large initial values of propeller pitch velocity end accelera- 
tion are indicated, followed by lesser values of propeller pitch velocity. 

4a 

dt 



Tire lasser values of 



are indicated in Figure 26. The initial values 



of id could not be measured. However, a plot of the change in 3 at 

this large initial ^ is given in Figiro 27. By inference, it may 
be theorized from Figure 27 that the lower the speed and the greater the 
amount of derivative control, the more extreme the demands or require- 
ments on the propeller control mechanism will be. This particular 
characteristic of response using derivative control should definitely 
bo investigated further. One v;ay measurements of this phenomena could 
be easily accomplished tJit h the differential analyzer is by changing 
the time base. 



It is realized that the turbo prop control problem is very 

\ 

complex and has not been solved to anyone's complete satisfaction to 
date, as is evidenced by the fact that the first production turbo urop 
system in the United States has only recently been installed in the 
Navy's XP5Y. This investigation has been based on possibly rudimentary 
considerations but a sincere attempt has been made to be as general as 
possible. It is realized that the investigation is not complete and 
that certain other approaches are desirable, even though not possible 
to investigate, due to lack of time. A discussion of some of these 
other considerati cns follows » 

A complete system analysis would necessarily consider the 
remaining variable not touched upon in this investigation, namely 
turbine bucket temperature. Such a system, as visualized by the authors, 
-would introdice this variable as a limit to fuel flow rate,, What the 
functional relationship between fuel flow rate and temperature should 
be or mitJit be has not been investigated. Hxwever, it is apparent that 
rate of change of fuel flow rate must be so limited, rogerdloss of 



throttle quadrant movement, to values such that temperature limita- 
tions will not be exceeded. In any event, the result would be less 
stringent than a step input of fuel flow, as i3 being used in this 
investigation, and results herein obtained should therefore be con- 
servative. 



governed exclusively by propeller blade angle. One problem envisioned 
In such an assumption would be in the case of simultaneous scheduling 
of RPM and fuel flow, which is in reality necessary as mentioned in 
the section on the derivation of uropeller data (III C). In this case, 
an increase in setting of fuel flow and RPM simultaneously would cause 
an initial decrease in blade angle due to the underspeei condition, 
resulting in an initially lower rather than higher torque and probably 
also in mare extreme oscillations. A discussion of this is ~iven in 
Ref. 2. 



only approach possible, although seemingly the most feasible and practi- 
cable. xn alternate approach as mentioned in Ref. (2), would cause fuel 
flow to bo a function of RPM only, with propeller blade angle controlling 
the remaining variables. 

Ho time lags, lost motion, or finite limitations on any of the 
variables were considered. This vrould be of importance in a physical 

system where, for instance, a •U were required instantaneously with 

d 't 

no time allowed for acceleration » but where the propeller governor 



problems would be caused by lost motion in linkages as in the propeller 



One basic assumption was to consider rotational speed 



Propeller blade fugle control of RRI is not, however, the 



could only provide a lesser 




at a finite acceleration. Other 




that the fuel pump response would be non-linear and asymptotic to both, 
original end new steady state values of fuel flow. 

Such phenomena would have the effect of non-linearising the 
system with attendant effect on response. It is considered that, once 
krDvm, those limitations could be easily set up on the differential 
analyzer. 

Another approach worthy of note would be to r. ake the dawning 
term in th e differential equation a function of ore or more variables, 
thereby makipg the dar.pirg variable. For instance, it might be desir- 
able to have low damping fqr large deviations from steady state and 
high danping far snail deviations from steady state. Such o system is 
described in Ref. 9, where the damping term is proportional to 



output angle • This results in the damping term approaching 

1/a at very small deviations from steady state »/ith lesser values at 
larger deviations, 'with such a system, very fast response can be ob- 
tained with little or no overshoot. Such a system could theoretically 
be set up for a turbo-prop system as hei'ein analyzed since at steady 




state 



k, Q - Ce w ^ ( Kf ) - ft»f. 



Also we know the t for this condition, >9 = . 

o i * 

Therefore, changing notation, « o = ^ (ut steady 



state). 



This means that tie toriue input to the system has a simple 



proportionality relat'ior. ship to a ficti ious , vfnich c-m then bo 




thereof could then be obtained, possibly by using a servo multiplier. 
This would e-pal the proportionality constant kg, getting back to the 
notation of this investigation. The complete damping tern would then 
be 



utilizing these principles would be difficult. Differentiation, inte- 
gration, addition, subtraction, multiplication, and division of variable 
would have to be accomplished in the control installation = Presentation 
of ti.is information, therefore, is intended to show merely what might be 
accorij li^ie d theoretically, using differential analyzer methods. 



which would be necessary to set up an actual control system on any of 
tne pre:rJ.ses investigated or discussed, since this investigation deals 
•with investigation, not with design. It is considered, however, that 
electronic differential analyzer non-linear components could be used in 
an actual installation, receiving measured variables like rotational 
speed, blade angle, torque, temperature, and fuel flow. These variables 
would be operated on, resulting in the formulation of correction signals 
As visualized by the authors , one such system 'would operate with blade 
angle controlling rotational speed and temperature limiting fuel flow. 
This /ould involve a differential analyzer component for the propeller 




It is not denied that the design of a workable system 



No effort has been made to indicate any physical hardware 



c ntrol and one for the fuel regulrtor. Torque would limit change in 



blade angle or fuel flow, or both, as deterrined to be desirable. 

» 

Since variable damping gives the best response, this might be incor- 
porated in the propeller control. A movement of the throttle control 
would then transmit new steady state values of KP.M end fuel flow to the 
propeller control analyzer component end the fuel regulator analyzer 
component, respectively. Response to these new steady state values 
would be accomplished in a minimum of time with no danger of exceeding 
any design limitations.? 

Such a systan would necessarily be complex and probably 
difficult to maintain and, at present, would be very difficult to design. 
Such a prospect should not be neglected, however, since it seems to have 



considerable merit 



SECTION V 



CONCLUSIONS i-2ID HECOi TJuNDATICNS 



1. The electronic differential analyzer presents an excellent means 
for the analysis of the requirements of e control system, such 

a3 that for a turbo-prop power plant, in which the complexity and 
interaction of tho variables makes mathematical analysis imprac- 
tical. 

2. The application of an electronic differential analyzer for the 
control of the turbo-prop power plant should be investigated. 

3. The natural-damping of the turbo-prop pov/er plant i3 insufficient 
for proper response without augmentation by artificial rneanSo 

4. By requiring the propeller pitch control to change pitch at a rate 
proportional to the engine offspeed md its integral, adequate 
system response can be obtained far any one operating condition 
by the proper choice of the prop pitch control proportionality 
cons tents. But, at other operating conditions, the response may 
not be satisf actoryo* 

5. By requiring the propeller pitch control to change pitch at a rate 
proportional to the engine offspeed, its integral, and its 



3P 



derivative, a consistent syston response at all operftinc 
conditions may be approached by te addition of the proper 
amount of derivative control. 

Oo The propeller pitch rate and accelor ti n should be limited to 
practical valuos to determine the effect on the response of the 
systems employed in this investig'ution. 

7. It is recommended that a further investigation of the control 
sys ten re uu i. re: art s i r c lud e : 

a. Turbine bucket temperature as a variable. 

b. The effects of control tine laps and dead time 
in Hie system. 

c. A study of other types of "rope Her oitch control 
to dotcrr.ine the cptinun t;qpe of Ofritiol, includ- 
ing non-linear requirements. 

d. « study of a systai in .hich a turbine fuel flou 
pcvurnor is used in place of the -propeller Hade 
an fie governor visualized in tl is investiy'- tion. 

e. A study of the control syslar response over vide 

% 

ranges of poser utilizing a scheduling procedure 
of engine soeed und fuel flow rate. 



*trrr 
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SECTION VII 









S/iMPLE CALCULATIONS 



Section III C. Selection of Propeller 

Determination of y«p anot Q Pe : 

D = prop, diameter = 1C ft. 
i = 11.35 



Assume: V = 300 IIFK = 440 ft/sec 

n = — =// 45 /?PA 1 = If. 

i 

£=0.002.378 SLUSs/ it i 

f 3 = 50° 

v A-wk^r z - 50Z 

Q = _ (O. e U5)(0-002378Xn*5 x X iol 

^ P p 2 TV n ~ z rr 



- 12,750 it. Lbs. 




12, 7SQ 
II. 35 



1 * er / 

/ /iec 



- l f )23 ft. Lbs. 



Section III E. Evaluation of Constants 

I = 129 lb. ft. 2 = 4.008 slug ft. 2 



I 

_E 

.2 

l 



70 



(llo35)‘ 



0.542 slug ft.‘ 



I * “§■ * 4.008 + 0.542 a 4.550 slug ft. 

e f 



Evaluation of Coefficients for Differential Equations 
Assuming: V q > 300 Uffl 



Jl & - 13,000 RPM 
vAT = WO lb/hr 



Operating Condition 



From I igure 4 : 



C 



. fiQe 
% Un< 



-.09 -ft. Lbs (bO 



w; 



From Figure 8: 



c - 3 -V 

c ^ e 



From Figure 11 



C =/^ e ) = 



re// . 

/nun. 


ITT, 


8: 




0.2 8 ft. lbs. . 


160 


RPM ' 


,27 T 


11: 




£4/ ft. Lbs. 


- i 



[zn) = ~ a 



U4- 



ft. Lbs. 



rtt */s 



Sec. 



= 2.(6 7 



11. Lbs. 









5 cJej \rees 



ft Lbs. 
deyr a 



For Stan 1 



= 0.6 

= 5 rad/sec 



Assuming: 




Ip 




4 -. 5 S sLuj 



a o = ur* a z = 5*{ +55) - //3. 6 



X 4 

siijy /r 
5ec l 



.j-2? /«»«,' = Z(0.t>) ) l(/l3.t)(+.tt)’ 



-- 2 . 7.3 



1 = £/? Ai 

O ' ,J ' 







y/3 .--^ = o. 8 r r 

nr./ 





-h 



k 2 * 






, Z7.S-2.UC7 +(- 0.7*4) 



rjL^, 



- /r e_ _ 



V 3 



= ft I if 7 



/ 2 £/ 



For Sten 2 



Assuming: £? = 0.6 

UJ ^ = 5 rad/sec 
5&).C. = 200/0 

a 2 [i + Z **-][ Z ‘+-fr] 

= [l -h Z ] £ 4.55-J = /3. 6 5" SLu 1 {t X 
= 5 Z 0*tS) = 



*1 



2 £ /^T 7 



s ,, c 

= y/. ? JL ±1ZIL. 

sec. 



a = Cp Hi 

o r ' 



o 

c i = -“^T 
v* 



34±£_ 

TxJT 



= Z.LL 



G 1 " ^*<-r7 

1 C/Ie 



e„ . + Cf p k *- 



ko = 



^ _ C Rjie + Ce^ e _ . 






/ zr.t 



> 



I* * $[** + %] 

[% P. cj fie * 7 «] 



k 3 ' 



. ^ 

U*>%)(+.ssl = ^ 07/ 









/2 8-J 



Section 17. Reduction of Data from Oscillograph Traces 



Example : 

Calculation of maximum torque for Step 2 with 
£ « 0<,6, - 5 rad/sec, 7 q = 300 MFH and 200 y D.C. : 

Maximum displacement in millimeters from trace => 13.3 MM 
Calibration = 20 ram/volt 
Attenuation = 100 

’Axltiplication factor = 1 (indicated on Figure 19) 

w * U3 - 3 ram) (volt/20 ram) (100) (1) = 66 0 5 volts 

= 66.5 ft. lbs. 



(See Table 7) 



SECTION VIII 



TABLES 

I. Perforr.ajice Data for the T. C. 100 
II. i reseller Data. 

III„ Differential Equation Coefficients 

IV. Data for Step 1 0 

V. Data for Step 2. 
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